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PACS 67 . 85 . -d - Ultracold gases, trapped gases 

PACS 03.65.Vf - Phases: geometric; dynamic or topological 

PACS 03 . 75 . Ss - Degenerate Fermi gases 

PACS 03 . 75 . Mn - Multicomponent condensates; spinor condensates 

Abstract. - We study Hubbard models for ultracold bosonic or fermionic atoms loaded into an 
optical lattice. The atoms carry a high spin F > 1/2, and interact on site via strong repulsive Van 
der Waals forces. Making convenient rearrangements of the interaction terms, and exploiting their 
symmetry properties, we derive low energy effective models with nearest-neighbor interactions, 
and their properties. We apply our method to F = 3/2, and 5/2 fermions on two-dimensional 
square lattice at quarter, and 1/6 fillings, respectively, and investigate mean-field equations for 
repulsive couplings. We find for F = 3/2 fermions that the plaquette state appearing in the 
highly symmetric SU(4) case does not require fine tuning, and is stable in an extended region of 
the phase diagram. This phase competes with an SU(2) flux state, that is always suppressed for 
repulsive interactions in absence of external magnetic field. The SU(2) flux state has, however, 
lower energy than the plaquette phase, and stabilizes in the presence of weak applied magnetic 
field. For F = 5/2 fermions a similar SU(2) plaquette phase is found to be the ground state 
without external magnetic field. 
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Ultracold atoms in optical lattices provide controllable 
quantum many body systems that allow to mimic con- 
densed matter (TJ [2] . They may in particular serve as 
quantum simulators of various Hubbard models [3], includ- 
ing those that do not have condensed matter analogues. 
Prominent examples include Hubbard models for bosons 
or fermions with high spin F. Experimental progress in 
studies of high F Bose-Eistein condensates [4] and Fermi 
gases (cf. (5 ) triggered a lot of interest in theoretical 
studies of such models. These studies go back to fun- 
damental questions of large N limit of SU(7V) Heisenberg- 
Hubbard model [6]; they have continued more recently in 
the context of ultracold atoms [7[|9]. These papers dis- 
cusses the interplay between the Neel, and valence bond 
solid (VBS), i.e. Peierls or plaquette ordering for antifer- 
romagnetic systems. Several other exotic phase should be 
possible of earth alkali atoms (cf. [To}{l3] ), where two or- 
bital SU(iV) magnetism, and even chiral spin liquid states 
were predicted. Several authors predicted also a variety 
of novel, exotic phases from effective (generalized Heisen- 
berg) spin Hamiltonians, obtained from spinor Hubbard 
models (cf. [14 15 ). A lot of effort was devoted to the in- 



vestigations of ID and ID ladder systems, where quantum 
effects are even stronger 16 . While for F = 1/2 Hubbard 
models quantum fluctuations suppress the conductor-Mott 
insulator transition [l7j, this is not the case for higher F, 
where dimer (Peierls) or valence bond crystal (Haldane) 
order, and in ladders even plaquette order are possible. 

Fermi systems with F = 3/2 were also intensively stud- 
ied [18]: first, because this is the simplest case beyond 
F = 1/2, second, because they can be realized with for in- 
stance with ultracold 132 Cs, 9 Be, 135 Ba, 137 Ba, and 201 Hg 
(for an excellent review see Ref. [l9]; Such systems ex- 
hibit a generic SO(5) or isomorphically, Sp(4) symmetry). 
In ID there exist the quartetting phase, a four-fermion 
counterpart of the Cooper pairing phase. In some situa- 
tions, counter intuitively quantum fluctuations in spin-3/2 
magnetic systems are even stronger than those in spin- 1/2 
systems. 

In this Letter we study the two-dimensional Hubbard 
model for F = 3/2 fermions with repulsive singlet and 
quintet interactions, and for F = 5/2 fermions in a 2D 
plane of its 3D parameter space. First, by rearranging the 
interaction terms and exploiting their symmetry proper- 
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ties, we derive low energy effective Hamiltonians. In con- 
trast to the standard approaches (see for instance fl4fl5 ), 
we do not use the spin representation, but rather keep the 
description in terms of fermionic operators. This allows 
us to formulate mean-field theory, somewhat analogous 
to slave-boson method [20], and show that the plaquette 
VBS state is stable in an extended region of the phase di- 
agram, in agreement with the predictions of Ref . |18[|19| . 
In the presence of weak applied magnetic field, however, 
the plaquette phase can be suppressed by an exotic SU(2) 
flux state. Moreover, for F = 5/2 fermionic atoms simi- 
lar SU(2) plaquette phase can be the ground state of the 
system without external magnetic field. 

Let us consider a system described by a Hamiltonian 
with nearest-neighbor hopping Hk in — — tJ3< 
and strong on-site repulsive interaction 
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(1) 



c \ a ( c i,a) are the usual creation (annihilation) operators 
of fermions with spin a at site i, and t is the hopping 
amplitude between the neighboring sites. Here, and in 
the following automatic summation over repeated Greek 
indices is assumed. The interactions depend on the spin 



of the scattering particles 21 : 
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where S is the total spin of the two scattering spin-.Fi 
and spin-F 2 particles. This means that the scattering 
processes can happen at different spin channels which are 
determined by the total spin of the scattering particles. 
Ps projects to the total spin-S subspace and gs is the 
coupling constant in the corresponding scattering chan- 
nel. Due to the on-site interaction the only contributing 
terms are either antisymmetric (as) or symmetric (s) for 
an exchange of the spin of the colliding particles depend- 
ing on their fermionic or bosonic nature. In the following 
we exploit this property of the on-site interaction which 
is also preserved for the effective strong repulsion model 
with nearest-neighbor interaction. 

Starting from the fundamental relation between the 
Ps projection operator and the product of the F 
spin operators for two spin-i* 1 fermions: (2FiF 2 )' = 

Y, s l S ( S + !) - 2F ( F + !)]' p s, the P s projector can be 
expressed as a degree of 2F polynomial of the F1F2 prod- 
uct: 

IF 

P S = Y, a SJ ( F i F 2)' 

1=0 



(3) 



for all S = 0, 1, . . . , 2F and with (FiF 2 ) = E. a s ,i are the 
coefficients of the expansion. Note that for a given value of 
S the Ps projector is either symmetric or antisymmetric 
in the spin indices of the scattering particles. This expan- 
sion is usually applied in order to express the high-spin 



two-particle interaction with effective multispin-exchange. 
In contrast, we will use the expansion of the projector op- 
erator in eqs. ([TJ and |2| in order to collect and treat ad- 
equately the two-particle interaction terms that describe 
different spin exchange and spin flip processes. In this case 
F denotes the three generators of the SU(2) Lie algebra in 
the appropriate representation: for high-spin fermions or 
bosons they are the SU(2) generators represented by the 
proper even/odd dimensional matrices. In case of pure bo- 
son or fermion system all processes take place only in the 
symmetric or antisymmetric part of the total spin space, 
respectively. Therefore, the following decomposition can 
be used: 



(FiF 2 )' = (FiF 2 ) 



(as) 



F1F2 
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and the symmetric and antisymmetric projectors can be 
expressed as follows: 



1=0 
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p ( s ] = E ^ (pW) 
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(5a) 
(5b) 



Here Nr as ) and iV( s ) denotes the number of antisym- 
metric and symmetric subspaces of the total spin space. 
The antisymmetric and symmetric part of an operator A 
can be constructed by the exchange of two spin indices: 

^ ( ° S C = Aa 4-<> and i Ai X,s = Aa J + A *n- lt 
is obvious that the above decomposition leads to the poly- 
nomials Eq. ([5]) having significantly smaller degree than 
Eq. ([3|. JV( as ) — 1 or iV( s ) — 1, respectively, determines the 
minimum degree of the polynomial of the product FiF 2 
which is equivalent to the interaction Eq. ([2]). 

Now let us apply the above procedure to a 2 dimensional 
F = 3/2 fermion system. In this case the interaction has 
to be antisymmetric therefore the only contributing terms 
are the total spin-0 (singlet) and the spin-2 (quintet) scat- 
terings: 
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At quarter filling there is only one particle on each site (i) 
on the average in general, and (ii) exactly if the on-site 
repulsion tends to infinity. Therefore for strong repulsion, 
the hopping can be considered as a perturbation, and the 
system can be described by an effective Hamiltonian with 
nearest-neighbor interaction. The effective model based 
on perturbation theory up to second (leading) order in 
the hopping t is the following: 



H eff = 



<i,j> 
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where V"f = J2i 
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gives the energy shift 



to the on-site energies due to the weak nearest-neighbor 
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hopping. Gs = —^t 2 /gs, (S = 0,2) is the new coupling 
constant in the spin-S* scattering channel. Since the effec- 
tive model preserves the symmetry of the on-site interac- 
tion, it remains antisymmetric for an exchange of two spin 
indices. Now the components of the F vector are the well 
known 4x4 spin matrices and Eq. (|5| has the following 
form: E( as ) = P a +P 2 , and (FiF 2 ) (as) = -15P /4-3P 2 /4. 
The effective Hamiltonian has the form 
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where a n = (5G 2 — Go)/4, and a s = (G 2 — Go)/3. 
The two-particle nearest-neighbor interaction terms are: 



Fig. 1: (Color online) The lattice was split into 4 sublattices 
(a) and due to this splitting the Brillouin zone shrank, (b) The 
shadow area depicts the reduced Brillouin zone. 

decoupled via a mean-field treatment by introducing the 
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J J „ „ r-p(as)i a ^ ar ,j /-p xp \(as) _ expectation values of the link operators (xi 1) and (Jj 

c i,a c j,P c j,S c tn L JLl J 7 ,o ' \ r l r 2) — j x] ._ /cf \ AT ____ j_i e'u tt.1;i 



c Ji7 [(FiF 2 ) (as) ]^,^, and their explicit spin de- 

<5a, 7 <5/3,o - 8 at g6prf, and 



pendence is: [E( ns )] a ' 

[(FiF 2 )(«)]^ = [FiL, 7 [F 2 ] w - [Fr] M [F 2 ]^ 7 . After 
straightforward calculations one arrives to the following 
form of the effective Hamiltonian: 



H eff = V {0) + E a n(n inj + xljXij - m) 



<hj> 



«,(s,s, • 
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where n,- 



and S, 



cL.F 



are the usual 



particle number and spin operators on site i, and 



Xi,i 



(10a) 
(10b) 



are introduced for the U(l), and SU(2) nearest-neighbor 
link operators, respectively. Note, that in general the 
SU(2) link operators do not satisfy the spin commutation 
relations, however, they clearly are related to the bond- 
centered spin. The competition between the spin and par- 
ticle fluctuations can be controlled by tuning of a n and 
a s . The effective Hamiltonian (|9| can be applied for less 
than quarter filled system too, provided the kinetic term 
is added to the Hamiltonian (|9|. contains the on-site 
energies and shifts the ground state energy only, so we do 
not consider its contribution. 

In the following we study the possible phases of the 
quarter filled system with the constraint J2 a c l a°i,a = 1 
— only single occupied sites are allowed due to the strong 
on-site repulsion. Due to this local constraint the Hamil- 
tonian is invariant under a rotation of the phase of the 
fermions at each sites. This means that the Lagrangian of 
the system £ = cj .9 T Cj i<T + H is invariant under the 
U(l) gauge transformation c;. CT — > Ci ya e^ i reflecting the 
local constraint for the particle number. 

Considering the Hamiltonian |9| the terms containing 
rii do not give contribution up to an irrelevant constant at 
quarter filling and the remaining 4-fermion terms can be 



and the spin operator (Si 
nian is: 

with 



Now the mean-field Hamilto- 



(11) 



H i,i = a «( (Xj,i) C L C J> + (Xij) c], a ^ ~ I (Xi,i) I 2 ) 

+ a s ({J jti ) c\ a ¥ a ^ Cj ^ + (Jij) ct ja F a)(9 Cj )( g - I (Ji,j) 



Si) cl F ai p Cj , p + (S 3 ) <i a F ai pCi,p - (S ; > (S 



Note that the mean-field Lagrangian also has to remain in- 
variant under the gauge transformation mentioned above. 
Thus, the link variables must transform as (Aij) —> 
(Ai tj ) e-^-^l 

The expectation values of the spin and link operators 
were determined self-consistently. Anticipating the ap- 
pearance of a plaquette phase similar to the ground state 
of the system for Go = G 2 , it is reasonable to split the lat- 
tice into 4 sublattices (see Fig. [T]) leading to the shrinking 
of the Brillouin zone to the quarter of its original size. We 
assume different values for the order parameters of the dif- 
ferent sublattices and of the alternating links as the only 
space-dependence of them. 

We have found the following gauge non-equivalent states 
with the following nonzero averages to be the solutions of 
the self-consistent equations: (a) Neel order: (Sj), (b) 
U(l) plaquette order: (Xi,j), an d (c) SU(2) plaquette and 
dimer order: (Sj), (Xi,j)i In the figures the nonzero 

expectation values are denoted by the following way: (Sj) : 
black arrow, (xi.j)'- black stripe, and (Jj t i): light blue 
stripe with arrow. The phase diagram of the system is 
shown in Fig. [2] If the effective interaction of the singlet 
channel is significantly stronger than that of the quintet 
channel, the dominant order is purely antiferromagnetic 
without any bond order. For a n < and a s > the spin 
and particle order compete with each other. In ref. 18 



a magnetically ordered dimer phase was suggested to ap- 
pear in this regime, however, we did not find any sim- 
ilar state to be the solution of the self-consistent equa- 
tions: when the antiferromagnetic order of the Neel phase 
is destroyed, plaquette order appears. The phase bound- 
ary is around Go ~ 1.9G 2 or equivalently — a n ~ 2.6a s . 
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(a) 



'' a n < 



a n = 



a s > 0, 



-a n S3 2.6a s 



;AFM 



U(l) plaquette 



a s = 



(b) 
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Fig. 2: (Color online) The phase diagram of the spin-3/2 
fermion system with strong on-site repulsion on 2D square lat- 
tice at quarter filling and the corresponding configurations: (a) 
Neel order (AFM), and (b) U(l) plaquette phase. 



In the U(l) plaquette phase the coupling constant a n al- 
ways dominates the interaction independently of the sign 
of a s . The nonzero U(l) links form boxes as shown in 
Fig. j2](b). One can define the U(l) plaquette variable as 
n = Xi,jXj,kXk,iXi,ii where i, j, k and I denote the sites 
of an elementary plaquette of the square lattice, and x 
is defined for nearest-neighbors only. The U(l) flux $ is 
defined by the phase of the plaquette. The plaquette vari- 
able and therefore the flux are invariant under the U(l) 
gauge transformations mentioned above. We have found 
two different gauge-non-equivalent states in the plaquette 
phase labeled by $ = and $ = it, respectively, and both 
states have the same energy. Note that our results are in 
good agreement with earlier results for the special SU(4) 
symmetric case (Go = G2), where similar box state was 
predicted with zero flux [8^ . We have found that this box 
state does not need fine tuning, it is the ground state in 
an extended, experimentally reachable region of the pa- 
rameter space. 

In the parameter region \Gq\ < 1 .9| G2 1 , and for a s 7^ 
we have found two other solutions on the top of the ground 
state having 10-15% higher energy: the SU(2) dimer phase 
and the SU(2) plaquette phase, where the latter corre- 
sponds to two gauge-non-equivalent states with different 
fluxes. Both the dimer and the plaquette phases have the 
same energy. In these states, in addition to weak ferro- 
magnetic order ((Si) < 3/2 and are equal for all of the 
4 sublattices) , both types of the link operators \ an d J 
have nonzero expectation values as shown in Fig. [3] (a), 
and (b). In the SU(2) plaquette phase the link operators 
with nonzero expectation values form plaquettes. These 
states are completely new RVS orders, therefore we make 
some notes about their basic properties and comment on 
their naming. Both states violate spin-rotation invariance, 
and the SU(2) dimer state — contrary to the SU(2) pla- 
quette phases — preserves the translational invariance by 
one lattice site in one spatial dimension. It is clear that J 
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SU(2) flux 
FM 



Fig. 3: (Color online) The magnetic field dependence of the en- 
ergy of the different states in units of t (the lines are only guides 
of the eyes), and the configurations in the SU(2) plaquette (a) 
and in the SU(2) dimer (b) phase. 



is not a member of SU(2) therefore it is reasonable to ask 
why do we use the terms SU(2) plaquette, flux or dimer for 
the states where the expectation value of J is nonzero? To 
answer this question let us consider the mean-field Hamil- 



tonian Eq. (11). The non-local part of the one-particle 



excitations appears in the Hamiltonian as 



( a n (xj,i) <W + a s (Jij) : 



a,/3 



a Cj,f>+H.c. (12) 



From this form it can be read that the excitations con- 
sist two branches with two different symmetries: (xj,i) 
relates to the U(l) excitations, while (Jj.j) F to the SU(2) 
excitations. In order to define the SU(2) flux let us in- 
troduce the new link parameter according to Eq. (12 1: 



Jij) F, with the usual inner product of the vectors 
in the 3 dimensional space of the generators F. Ui.j is a 
member of SU(2) and a 4 x 4 matrix for F = 3/2 fermionic 
atoms and the same holds for the SU(2) plaquette variable 
n sc/(2) = U itj Uj, k U kt iUi ti . The flux * passing through the 
plaquette defined by the form: Tl su ( 2 } = e** F . In order to 
determine the ground state it is worth to express the mean- 
field Hamiltonian with the Ji.j operators, while the excita- 
tions and the SU(2) flux can be expressed with Ui.j. Note 
that the SU(2) plaquette H su ^ is also invariant under the 

r . e i<t>i 

Consid- 



U(l) gauge transformation defined above: Cj CT — 

(Xi,j) -> (Xi,j) e^"H and U id -> U^e*^^ 



ering the definition of Uij, the last relation is obviously 
equivalent to the transformation J,j — > J^ j-e^w - ^'). 

The SU(2) phases can patently claim to great interest, 
but they are suppressed by the U(l) plaquette state. Nev- 
ertheless, since the SU(2) flux, as well as the SU(2) dimer 
order coexist with ferromagnetic order, it can be expected 
that weak magnetic field does not destroy the SU(2) order, 
but it can stabilize that. To check this let us investigate 
the energy of these states in the presence of external mag- 
netic field h taken into account as a Zeeman term in the 
Hamiltonian: 



H" 



H 



MF 



(13) 
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Note that for strong magnetic field the quadratic Zeeman 
term can become important. At this point let us suppose 
that this term can be neglected and we will check the 
validity of this assumption at the end of the calculations. 

The magnetic field dependence of the energy of the 
SU(2) flux state compared to the U(l) plaquette state and 
to the ferromagnetic state is shown in Fig. [3] for a typi- 
cal value of the couplings in units of the nearest-neighbor 
hopping. The U(l) plaquette phase remains the ground 
state in the presence of nonzero, but very small magnetic 
fields. The SU(2) plaquette state, as well as SU(2) dimer 
order, have the lowest energy for higher value of the ap- 
plied magnetic field h, and are the ground state of the sys- 
tem in an extended region of the phase diagram. Both the 
SU(2) dimer and the plaquette states either or ir flux 
have the same energy. In the presence of even stronger 
magnetic fields ferromagnetic order suppresses any other 
order in the system. Now, let us check the validity of our 
assumption of neglecting the quadratic Zeeman coupling. 
The linear Zeeman energy (in h unit) is given by the Lam- 
our frequency: lj^ — gFfJ-BB, and the quadratic Zeeman 
energy is: u q — jj^j. Here u>hf is the hyperfine split- 
ting energy, gp is the gyromagnetic factor and /is is the 
Bohr magneton. If ^ 3> 1 or equivalently ^> 1, the 
quadratic Zeeman effect can be neglected. We measure 
the magnetic field in units of the hopping parameter t. In 



-G 



optical lattices t — u)r-£, e ^' , where £ = (Vq/Vr 



a/4 



V 



is the potential depth, and ujr is the recoil energy, t has 
a maximum at Vo rs cjr, where t ~ ujr. lor is typically in 
the order of 1 - 100 kHz (however, wr/2tt - 400.98 kHz 
for 9 Be). We have found that the SU(2) flux state has 
the lowest energy if the linear Zeeman energy is around 
u>l ~ 0.1 — It which means 0.1-100 kHz. The hyperfine 
frequency is in the order of 1-10 GHz, so uihf/^L is in the 
order of 10 4 — 10 7 . Therefore, the quadratic Zeeman effect 
can be neglected compared to the linear Zeeman term for 
magnetic fields which are sufficient to stabilize the SU(2) 
phases. There are some atoms which could be promising 
candidates for realizing experimentally spin-3/2 fermion 
systems and have much less recoil energy than the above 
mentioned 6 Be, namely 1 — 10 kHz which in our case corre- 
sponds to very small magnetic field: 10~ 5 — 10~ 4 G. With 
these atoms the experimental realization of the SU(2) or- 
ders demands quite strong magnetic shielding. 

Similar analysis can be easily made for F — 5/2 
fermions at 1/6 filling for the special values of the cou- 
pling constants G4 ~ (— 7 Go + 10G2)/3. G4 is the cou- 
pling of the interaction with 9-fold spin multiplicity that 
appears in the Hamiltonian Eq. (|6| for spin-5/2 system 
in addition to the singlet (Go) and quintet (G2) scat- 
terings. In the plane of the parameter space defined by 
G4 = (— 7Go + IOG2V3, the structure of the Hamiltonian 
is exactly the same as Eq. (Ill, there is no term con- 
taining higher order of the product F1F2, and the cou- 
plings take the values: a n = (— 23Go + 35G2)/12 and 
a s = (— Gq + G2V3. Note, that now the components of 



a n > 0/ «s > 0,-' 

/an < /' 
an = / Sa 3 < 
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AFM 



a s = 
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SU(2) plaquette 
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Fig. 4: (Color online) The phase diagram of the F = 5/2 
fermion system with strong on-site repulsion on 2D square lat- 
tice at 1/6 filling and the configurations in the U(l) quasi- 
plaquette phase (a), and in the SU(2) plaquette phase (b). 



the F vectors are the 3 generators of the SU(2) Lie alge- 
bra in 6 x 6 representation. During the calculations we 
have used the same 4-sublattice ansatz as in case of the 
F = 3/2 fermions, because the more suitable 6-sublattice 
ansatz does not respect the symmetries of the original 
model on square lattice. Similarly we have determined 
the solutions of the self-consistent equations for the (S j) , 
(Xi,j) & n d (Jj,j) expectation values and the corresponding 
energies. The ground states of this system are shown in 
Fig. [4] In case of dominant singlet scatterings the ground 
state is purely antiferromagnetic, at least while a n > 0. 
A weak negative a n seems to lead to an instability in the 
system, but we could not find any stable solution of the 
self-consistent equations in this narrow region. Further 
decreasing |Go|, a quasi-plaquette phase appears. In this 
phase the expectation value of the U(l) link operator (Xi,j) 
is non-zero everywhere, but stronger and weaker links al- 
ternate forming a weak plaquette structure. The flux pass- 
ing through the plaquette is zero and there is no spin order 
in this phase. For even weaker singlet coupling (increas- 
ing the value of \a s \) weak ferromagnetic order appears 
in addition to the plaquette order (the SU(2) plaquette 
phase in Fig.|4|. Here the plaquettes are formed not only 
by the alternating zero and non-zero U(l) link operators 
(Xi,j), but the SU(2) operators (Jjj), too. The flux pass- 
ing through the plaquettes remains zero. This means that 
while for F = 3/2 fermions the SU(2) plaquette phase 
can be the ground state of the system only by applying 
external magnetic filed, for F = 5/2 fermions a similar 
SU(2) plaquette order (with zero flux) is the ground state 
of the system for weak singlet couplings. It is difficult to 
determine precisely the phase border between the quasi- 
plaquette and the SU(2) plaquette phases. This is due 
to the fact that the two orders start to compete around 
2Go ~ G2: both are stable solution of the self-consistent 
equations but approximately with the same energy. The 
energy difference between the two phases increases slowly 
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for decreasing |Go|, and around 3Go ~ G2 reaches the 
4 — 5%. Note, that in the same parameter regime we have 
found another stable solution: a plaquette phase with it 
flux and with stronger ferromagnetic order, however (Sj) 
remains smaller than 5/2. In this state the SU(2) order 
parameter dominates the link variables: the value of (Jij) 
is twice than that of the corresponding (Xi,j)- The energy 
of this spin ordered SU(2) plaquette state with n flux is 
higher by about 5% than the one with zero flux. 

Finally let us discuss the validity of our mean-field re- 
sults. In two dimensions the mean-field solutions are ex- 
pected to provide qualitatively reliable results. However, 
fluctuations around the mean-field results can become rel- 
evant, especially in the following two cases: when the sys- 
tem is close to the phase boundary of a continuous tran- 
sition, the correlation length becomes very large and fluc- 
tuations cannot be neglected. The present work focuses 
on the qualitative description of possible exotic states of 
matter. The states discussed above are the ground state of 
the system in an extended region of the parameter space 
and far enough from the phase boundaries expectedly they 
can be realized experimentally. Fluctuations could also 
be important if the mean-field solution is degenerate and 
the states in the degenerate subspace are not separated 
by energy barrier. In our case e.g. the SU(2) dimer and 
plaquette states have the same energy but they are sepa- 
rated by energy barrier: one cannot arrive from one to the 
other with continuous (link by link) deformation without 
increasing the energy. These types of mean-field solutions 
are expectedly not effected by small fluctuations. 

To summarize, we have used a decomposition of the to- 
tal spin space into its symmetric and antisymmetric part 
with respect to the exchange of two spin indices of the high 
spin scattering particles. This decomposition was used for 
strongly repulsive systems to derive effective low energy 
Hamiltonians. This task was achieved remaining within 
the two-particle representation. The main advantage of 
the treatment is that it does not require to introduce com- 
plicated effective multiparticle/multispin interactions, but 
relies only on rearrangements of the usual two-particle in- 
teractions. The effectiveness of the treatment does not 
depend on the statistics of the considered particles, and 
it allows us to identify the different processes in the spin 
channel within the concept of site and bond spin. Apply- 
ing this method to F = 3/2 fermions, we have determined 
the ground state phase diagram of the system on mean- 
field level to complete the earlier results known for some 
regimes of the couplings. We have found that the VBS 
state is the ground state in an extended region of the phase 
diagram, while in the presence of weak magnetic field an 
exotic SU(2) flux state has the lowest energy. We have also 
made some similar calculations for F — 5/2 fermions in the 
plane determined by the condition G4 = (— 7Go + 10G2)/3 
in the 3-dimensional parameter space of the coupling con- 
stants. We have found that novel, exotic SU(2) plaquette 
phase, similar to one predicted for spin-3/2 fermion sys- 
tem in the presence of external magnetic field, can be the 



ground state of the system with zero flux, without applied 
magnetic field. 
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